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O Abstract 

^-j- We prove that the entanglement created in the low-energy scattering of two particles in two dimensions is given by 
a universal coefficient that is independent of the interaction potential. This is strikingly different from the three 
dimensional case, where it is proportional to the total scattering cross section. Before the collision the state is a 

Oh 

{ product of two normalized Gaussians. We take the purity as the measure of the entanglement after the scattering. We 
give a rigorous computation, with error bound, of the leading order of the purity at low-energy. For a large class of 
potentials, that are not assumed to be spherically symmetric, we prove that the low-energy behaviour of the purity, 
V, is universal. It is given by V = 1 — j^^jj^i £, where o is the variance of the Gaussians and the entanglement 
coefficient, £, depends only on the masses of the particles and not on the interaction potential. The entanglement 
depends strongly in the difference of the masses. It takes its minimum when the masses are equal, and it increases 
rapidly with the difference of the masses. 



1 Introduction 

In this paper we consider the low-energy scattering of two particles without spin in two dimensions. The interaction 
between the particles is given by a general potential that is not required to be spherically symmetric. Before the 
scattering the particles are in an incoming asymptotic state that is a product of two s. After the scattering the 
particles are in an outgoing asymptotic state that is not a product state. The problem that we solve is to compute 
the loss of purity of one of the particles, due to the entanglement with the other, that is produced by the collision. 

In the configuration representation the Hilbert space of states for the two particles is V. :— L 2 (R 6 ). The dynamics 
of the particles is given by the Schrodinger equation, 

d 

ih—<fi(x 1 ,x 2 )=H(p(x 1 .x 2 ), (1.1) 



*PACS Classification (2010): 03.67.Bg, 03.65.Nk, 03.65. Db. Mathematics Subject Classification(2010):81P40, 81U05, 35P25. 
tResearch partially supported by CONACYT under Project CB-2008-01-99100. 
t Fellow, Sistcma Nacional de Investigadores. 



f 



where the Hamiltonian is the following operator, 

H = ff + y(x 1 -x 2 ), (1.2) 

where Hq is the free Hamiltonian, 

By h it is denoted Planck's constant, rrij,j = 1,2, are, respectively, the mass of particle one and two, and Ay, is the 
Laplacian in the coordinates Xj,j = 1,2, of particle one and two. The potential of interaction is multiplication by 
a real- valued function, V(x), defined for x G K 3 . We suppose that the interaction depends on the difference of the 
coordinates xi — x 2 , but we do not require the spherical symmetry of the potential. We assume that V satisfies mild 
assumptions on its regularity and its decay at infinity, namely that. 

ASSUMPTION 1.1. 

(1 + |x|)^y(x) G L 2 (R 2 ), for some /3 > 11. (1.4) 
Under this assumption if is a selfadjoint operator. 

We also suppose that at zero energy there is neither an eigenvalue nor a resonance (half-bound state), for the 
Hamiltonian for the relative motion H re i := — + V(x) where x e R 2 is the relative distance and m is the reduced 
mass. A zero energy resonance (half-bound state) is a bounded solution to H re iip = that is not in L 2 (M. 2 ). See [T] 
for a precise definition. For generic potentials V there is neither a resonance nor an eigenvalue at zero for H re i. That 
is to say, if we consider the potential XV with a coupling constant A, zero can be a resonance and/or an eigenvalue for 
at most a finite or denumerable set of A's without any finite accumulation point. 

We study our problem in the center-of-mass frame. We consider an incoming asymptotic state that is a product 
of two normalized Gaussian, given in the momentum representation by, 

¥V Po (pi,p 2 ) :=p Po (pi)¥>-p (p 2 ), (1.5) 

with, 

^(Pi)-^^-^) 2 /- 2 , (1.6) 
where pj, i = 1,2 are, respectively, the momentum of particles one and two. 



In our incoming asymptotic state ( 1.5 1 particle one has mean momentum p and particle two has mean momentum 
— Po- Both particles have the same variance, a, of the momentum distribution. As we suppose that the scattering 
takes place at the origin at time zero, the average position of both particles is zero in the incoming asymptotic state 



(1.5 1. The outgoing asymptotic state of the two particles, v?out,p : after the scattering process, is given by, 



¥>out,p (Pi:P2) := (S(p 2 / 2 w)^i„, Po ) (Pl,P2)- (1-7) 
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Here, p := — T " 2 pi ^ — p 2 is the relative momentum, m := m,\ m?,/ {m\+m<i) is the reduced mass, and S(p 2 /2m) 



is the scattering matrix for the relative motion. Observe that in the state (1.5) the mean relative momentum of the 
particles is equal to p . 

The measure of entanglement of a pure-bipartite state that we use is the purity of one of the particles. Namely, 
the trace of the square of the reduced density matrix of one of the particles, that is obtained by taking the trace on 
the other particle of the density matrix of the pure state. Note that the purity of a product state is one. 

The purity of tp out , Vo is given by, 

V{ipont,p ) = / rfPldpirfP2rfP2'^out,p Q (Pl,P2)'^out,po(pi,P2)'^out,p Q (Pi,P2)'^out,po(Pl,P2)- (1-8) 

JM. 12 

Since the relative momentum, p, depends on pi and on p 2 , <p O ut,p is no longer a product state and it has purity 
smaller than one. This implies that the collision has created entanglement between the two particles. 

To be in the low-energy regime the following two conditions have to be satisfied. 1. The mean relative momentum 
Po has to be small. 2. The variance a has to be small. Note that if a is large the incoming asymptotic state </5i n ,p 
has a big probability of having large momentum, even if the mean relative momentum p is small. 

Let us designate by <p- m the incoming asymptotic state with mean relative momentum p = 0. The corresponding 
outgoing asymptotic state is (p ou t ■— S(p 2 /2m)(p- m . 



We denote by 



777 " 

m:= ,1=1,2, (1.9) 

mi + m 2 



the fraction of the mass of the i particle to the total mass. 



In Theorems |3.3| and |3.6| in Section 3 we give a rigorous proof of the following results on the leading order of the 
purity at low energy. 



7% out , P0 ) = V M + J^l O ( lHa/h l lpolM2 ) , asa/ft + |po|/» - 0, (1.10) 

V ^ = 1 - (m^jw 8 ^ + i\wiw) ' ™ a/h -> °' (L11) 



where £(ni) is the entanglemement coefficient 

2ir 2 
1 + (2 Ml - 1) 



g 0*i) : = 1 ■ ^ iv» I 1 + V 71 + Wx - !) 2 1 - 2 Mi) + J(l- Ml, Mi)], (1-12) 

J 7T 



with 

^(Mi,M2) := J 'dq 2 [ J dqi Exp[-i(/i? + /i|)(q x + q 2 ) 2 - (m2Qi - MiQ2) 2 - qf/2] 

VlMi ~ M2I |qi + qal iMaQi - MiQ2|)] 2 , 
where Iq is the modified Bessel function [2J. 



(1.13) 
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Observe that S(fJ-i) = S(l — fix), as it should be, because V(<p ou t) is invariant under the exchange of particles one 
and two. 

Note that J(l/2, 1/2) = tt 3 . In the appendix we prove explicitly that J(l, 0) = 16.6377. For [i x e [0, 1] \ {1/2, 1} 
we compute J (pi, 1 — Hi) numerically using Gaussian quadratures. 

The entanglement coefficient E (/ii) is universal in the sense that is independent of the interaction potential. Of 
course, this is only true if the potential is not identically zero, see the low-energy estimate for the scattering matrix 



given in (2.11). This is strikingly different from the three dimensional case that we previously studied in [3J, where 
the entanglement created by the collision is proportional to the total scattering cross section. 

Table 1 and Figure 1 show that -as in three dimensions [3J- the entanglement coefficient depends strongly in the 
difference of the masses. It takes its minimum for fii = 0.5, when the masses are equal, and it increases rapidly 
with the difference of the masses, as fix tends to one. This implies that in experimental devices intended to create 
entanglement by collisions it is much more convenient to use particles with a large mass difference. There is a physically 
intuitive reason for this: in the collision of a heavy particle with a light one, the trajectory of the light particle is 
strongly changed, and there is large exchange of quantum information between the particles, what produces a large 
entanglement creation. 

It is also interesting to remark that in the scattering of a particle with a large mass and a particle with a small 
mass we can assume that the trajectory of the large particle is not affected by the interaction, i.e. that, to a good 
approximation, it follows a free trajectory, and that the small particle feels a (external) interaction potential centered 
in the position of the large particle. However, the trajectory of the small particle will be strongly affected by the 
interaction, what will produce exchange of information between the particles, leading to the creation of entanglement 
between them. To evaluate this entanglement it is, however, necessary to take into account the degrees of freedom of 
both particles, as we do to compute the purity. 

Besides its intrinsic interest, there are many reasons why it is important to study the creation of entanglement 
in scattering processes. From the conceptual point of view scattering is probably the simplest way to entangle two 
particles. Before the scattering, in the incoming state, the two particles are in a pure product state where they are 
uncorrelated. As they approach each other they become entangled by sharing quantum information between them. 
After the scattering, they are far apart from each other, but they remain entangled in the outgoing asymptotic state, 
that is not anymore a product state. Scattering is a basic dynamical process that is essential across all areas of physics. 
Furthermore, scattering interactions are fundamental at all scales and there is a large variety of scattering systems. As 
is well known, entanglement is a central notion of modern quantum theory. It is the fundamental resource for quantum 
information theory and quantum computation. It is a measure for quantum correlations between subsystems. In the 
case of bipartite systems in pure states, entanglement is a measure of how far away from being a product state a 
pure state of the bipartite system is. Currently it is well understood that entanglement in a pure bipartite quantum 
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state is equivalent to the degree of mixedness of each subsystem. For information about this issue see [5], [BJ, [7J. 
Moreover, the study of entanglement creation in scattering is interesting for a many other reasons. For example, for 
the implementation of quantum information processes in physical systems where scattering is central to the dynamics, 
like ultracold atoms and solid state devices. Moreover, the study of entanglement in the scattering of particles requires 
quantum information theory with continuous variables and mixed continuous-discrete variables. See [7J for a review 
of this topic. It is possible that scattering will provide a new perspective to quantum information theory. Finally, 
entanglement creation is important to the theory of scattering itself, because it poses new problems that can shed new 
light and new points of view in the study of scattering processes. 

For previous results in the generation of entanglement in scattering processes in one dimension, mainly for potentials 
with explicit solution, see [4], [8], and the references quoted there. For the three dimensional case see [3]. Furthermore, 
[2], [lO], [H], and the references quoted there, consider a system of heavy and light particles. They study the asymptotic 
dynamics and the decoherence produced on the heavy particles by the scattering with light particles in the limit of 
small mass ratio. This problem is different from the one that we discuss here. The loss of quantum coherence induced 
on heavy particles by the interaction with light ones has attracted much interest. See for example [12) . and |13j . 

The paper is organized as follows. In Section 2 we study the low-energy asymptotics of the scattering matrix for 
the relative motion of the particles. In Section 3 we prove our results in the creation of entanglement. In Section 4 we 
give our conclusions. In the Appendix we explicitly evaluate integrals that need. Along the paper we denote by C a 
generic positive constant that does not necessarily have the same value in different appearances. 

2 Scattering at Low-Energy in Two Dimensions 

We denote by Ti. := L 2 (M. ) the state space in the momentum representation. The momentum of the particles one and 
two are, respectively, pi, P2. It convenient to take as coordinates in the momentum representation the momentum of 
the center of mass and the relative momentum, 

Pern := Pi + P2, (0 n 

,_ ro 2 pi-miP2 

" ' 1111+1112 

The state space in the momentum representation factorizes as a tensor product, 

H = H cm ®H ieh (2.2) 

where H cm — L 2 (R 2 ),H rc \ := L 2 (R 2 ) are, respectively, the state spaces in the momentum representation for the 
center-of-mass motion and the relative motion. 

Since the potential depends on the difference of the coordinates of particles one and two the scattering matrix for 
the system decomposes as the tensor product I crn C3> S(p 2 /2m) of the identity on "H cm times the scattering matrix for 
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the relative motion, <S(p 2 /2m), in 7^ re i, where m is the relative mass, 

m := (2.3) 
mi + m 2 

The scattering matrix <S(p 2 /2to) is a unitary operator in L 2 (S X ) for each p 2 /2m G (0, oo), where we denote by S 1 
the unit circle in M 2 . 

We introduce some notation that we need. We denote v :— \/\V(x)\ . Let P, Q be the projector operators in 
L 2 (R 2 ), 

P:= -«(x) (•,«), Q:=l-P, (2.4) 
a 

where, 

a := / |V(x)|dx. (2.5) 

Furthermore, 

= f i, ify(x)>o, 

■ | _j if y( x ) <0 . ^- Dj 
By Moo we denote the integral operator with kernel, 

M 00 (x,y):=C/(x)5(x-y)-i-^ U (x)ln^^=^) «(y), (2.7) 

where 7 is Euler's constant. Moreover, by Noo we denote the integral operator with kernel, 

1 9m 

7V ,o(x, y) := C/(x) <5(x - y) - «(x) In (|x - y|) „(y), (2.8) 

and 

A, := (QMqoQ)" 1 , a bounded operator, QL 2 (R 2 ) -> QL 2 (K 2 ). (2.9) 

The assumption that is neither a resonance nor an eigenvalue for H re i precisely means that (QM00Q) is invcrtiblc 
on QL 2 (R 2 ) with bounded inverse. 



Finally, we designate, 



Y Q (y) := (2.10) 

V 27T 



For X, Y Banach spaces we denote by B (X; Y) the Banach space of all bounded linear operators from X, into Y. 
In the case X = Y we use the notation B(X). By TrA we designate the trace of the operator A. 

THEOREM 2.1. Suposse that Assumption \ is satisfied and that at zero H re i has neither a resonance (half-bound 
state) nor an eigenvalue. Then, in the norm of B (L 2 (S 1 )) we have for \p/h\ — »• the expansion, 

S( P 2 /2 m) = / + «r ^ E + (Wn 2 - 7 + I) - £ ) ^-^^ E + O (j^JJ^) , (2.11) 



G 



where I is the identity operator on L 2 (§ 1 ), 



£:= (;Y )Y , 



and a is the scattering length defined by 

1 2tt 



Tr [PN 00 P - PMooQ^ooMoo + PM 00 Q] . 



(2-12) 



(2-13) 



a a 

Proof: Let us denote by <Si(A) the scattering matrix for the Hamiltonian Hi := —A + %«rV(x). It follows from an 



elementary argument that, 



Furthermore [14] . 



where T(A) is the trace operator, 



S(p 2 /2m)=S 1 ((p/H) 2 ). 



S 1 (\)=I-2iri 2 ^T(\)v (M^y 1 «r*(A), 



(r(A) V ) (*/) := 4? ^ / e- iyXl "Xx)dx, .eS 1 , 

V2 27T J R 2 



and 



M(A) := ff7 + |J«(-A - A - iO^ 1 ^ 



has a bounded inverse in i 2 (]R 1 ). 
Moreover, for all /ei^S 1 ), 



,(x) |P-is] I 



1 



L 2 (M 2 ) 



Furthermore, by Schwarz's inequality, 

Hx)|x|||i 2(R2) = ||nx)|x| 2 ||, !; ||( I +|/|) T(x) 

and it follows that, 



|(i + ixirHi L2(R2) <c, 



u(x)T*(A) = -^u(x)E + 0(\/A), A^O, 
v2 



in the operator norm in B(L 2 (S ), L 2 (R 2 )). Taking the adjoint in both sides of (2.19) we obtain that, 



1 1 



TOur^) = ——{.,v) + 0(V\), A^O, 



in the operator norm in B{L 2 (R 2 ), L 2 (S X )). By (2.15), (2.19), (2.201 



in the norm of B (L 2 ^ 1 )). Equation (|2.11[) follows from (|2.21[) and Theorem 6.2 of [I]. 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



□ 



The low-energy expansion (2.11) was previously proved by [15, in the case of exponentially decreasing potentials 
such that J V(x) ^ 0. 
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3 The Creation of Entanglement at Low-Energy 



In what follows we consider a pure state of the two-particle system. The wave function in the momentum representation 
is given by </?(pi,p 2 )- We designate by p(ip) the one-particle reduced density matrix with integral kernel, 



p(</>)(Pi,Pi) ■= J ¥>(Pi,P2Mp'i,P2) dp 2 - 



The purity, V((f), is given by, 

V(<p) :=Tr(p 2 ) = J dpidpidp 2 dp2^(pi,p 2 )¥'(pi,P2)^(pi,P2)¥'(Pl>P2)- (3- 1 ) 

As is well known [7J [THl [17] , the purity is a measure of entanglement that is closely related to the Renyi entropy of 
order 2,— In Tr(p 2 ). Moreover, It is trivially related to the linear entropy, Si, as Sl = 1 — V. Clearly, It satisfies 
< V < 1 if tp is normalized to one. Furthermore, it is equal to one for a product state, ip — Pi(pi) <^ 2 (p 2 )- As we 
will show, the purity can be directly computed in terms of the scattering matrix. For this reason it is a measure of 
entanglement that is convenient for the study of entanglement creation in scattering processes. 

We consider an incoming asymptotic state, in the center-of-mass frame, that is a product of two normalized 
Gaussian wave functions, 

¥V,p (Pi>P2) := < y5 Po (Pi)^-po(P2) J (3-2) 

where 

^po(Pi)-p^r72^ (pi - po)2/2CT2 - (3-3) 



Note that in our incoming asymptotic state (3.2) particle one has mean momentum p and particle two has mean 
momentum — p . The variance of the momentum distribution of both particles is a. We assume that the scattering 
takes place at the origin at time zero, and for this reason the average position of both particles is zero in the incoming 



asymptotic state (3.2). Observe that by (2.1) the mean value of the relative momentum in the state (3.2) is equal to 
Po- 

Since the incoming asymptotic state Pi n ,p is a product state its purity is one, 

P(Wn,Po) = 1- (3-4) 

The outgoing asymptotic state of the two particles, (p ou t,p Q -after the scattering process is over- is given by 

Vout, Po (pi,p 2 ) := (5(p 2 /2m)^i„ : p ) (pi,p 2 ). (3.5) 

As the relative momentum p depends on p! and on p 2 , p ou t,p a is not a product state, and then it has purity smaller 
than one. This implies that the scattering process has created entanglement between the two particles. 

Below we give a rigorous computation of the leading order of the purity of the outgoing asymptotic state f ou t tPo 
in the low-energy limit for the relative motion. To be in the low-energy regime the following two conditions have to 
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be satisfied. 1. The mean relative momentum p has to be small. 2. The variance a has to be small. Note that if a is 
large the incoming asymptotic state ^in.po bas a big probability of having large momentum, even if the mean relative 
momentum p is small. 

Let us introduce some notations that we use later. 

We denote by tp- ul the incoming asymptotic state with mean value of the relative momentum zero, 



where, 



¥>in(Pi,P2) := tp(pi) ip(p 2 ), 



¥>(P) 



( C r 2 7T) 1 /2 

and by tp out the outgoing asymptotic state with incoming asymptotic state ip- m 

Pout(Pi,P2) := (S(p 2 /2m)cp ia ) (pi,p 2 ). 



(3.6) 



(3.7) 



(3.8) 



Recall that, 



Mi 



m l + 1 Jl 2 

is the ratio of the mass of the i particle to the total mass. 



,* = 1,2, 



It follows from (2.1) that, 



Pi = MiPc 



P2 = A*2Pcro - P 



We prepare some results that we will use. It follows from (3.2 3.3| 3.10| 3.11 1 that 

, n . - _ p-(M?+M2)Pcm/ 2<T2 P -(P-PO) 2 /0- 2 P -(A"1-M2)pcm-(P-P0) 

G lr K 

REMARK 3.1. For some positive constant 8, 



bin I < ~2~ e 



-<5(pL+P 2 )/2^ 2 



Proof: Note that [i\ + [i\ = \ (l + (/ii - /i 2 ) 2 ). Then, for a > 0, 

(Ml + A^Pcm/ 2 ^ + P 2 /(T 2 + (/Ul - M2)Pcm ' P > 

(i + (Ml - M2) 2 (| - f )) ^ + S (1 - sk) > %L + P 2 )/2a 2 , 



(3.9) 



(3.10) 
(3.11) 

(3.12) 
(3.13) 



(3.14) 



provided that we choose a so that, < 5/2 < min[| + (/Xi — /i 2 ) 2 (| — f ), (1 — -p)]- The remark follows from (3.12) 



and (3.14). 



PROPOSITION 3.2. For any a, /3 > 0, 

(ln(2 + ]p|A)r 



(l + \ln(\p\/hW ^ n ' P ° 



\n(a/h+\p \/h)f 



&sa/h + \po\/H—> 0, 



(3.15) 
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(ln(2+|p|A)) a 



I MM A)!" 

Furthermore, uniformly for |po|/<7 in bounded sets, 
(ln(2 + |p|/ft)) 



O 



\M<r/K)\ p J' 



as a /H—> 0. 



(i + IMIpI/s)!)" 



(^in,p - ¥>in) 



|lm>/ft+|p |/ft)|V 



asa/H + |po|/ft — > 0. 



Proof: By (3.13) we have that 



(ln(2+|p|/fi)) Q 
(l + |ln(|p|/ft)|)^ 



where, for some 1 > 7 > 0, 



h ■= 



(ln(2+|p + p |/fi))° V 



,-<5(pL+p 2 )/* 2 



'|p|/ CT >l/(<r/ft+|p |/S)T \(1 + |ln(|p + Po|A)|)' 3 / CT 4 7T 2 

7V = 1,2,- ••, and 



dPcm rfp < C N (a/h + I po| /ft) 



J a ■= / / (ln(2+|p + po|/fi))° \ 2 e _ 5(pL+p2)/ ^ ^ dp 

/|p|/<r<l/(<r/ft+|pa|/ft)-r V i 1 + I ln ( IP + Po I ) /3 / ^TT 2 

Moreover, for |p|/<r < l/(cr/ft+ IpoIA) 7 and (V/ft + Ipol/fr) 1 " 7 < 1/2, 



< 



(l + |ln(|p + p |/ft|)| " (l + (l-7)|ln(Wfi+|p |/a) 2Vd-7)|' 



and then, 



h = ° 1 UnW^PoWFj ' '/*+|Pol/*-0. 



Equation (3.15) follows from (3.191 and (3.21) 



In the same way, it follows from (3.13) that 



(ln(2+|p|/ft))« 



|ln(|p|/ft)P 



< h + h, 



where, for some 1 > 7 > 0, 



h 



(ln(2 + |p|/ft))<*V 1 



|p|/<r>l/(a/ft)-r V I MIpI /K)\ J ^ 



-«(pL+p 2 )/^ 2 



dp cm dp<C N (a/h)", N = 1,2, 



and 



'| P |/<7<i/(cr/fc)-» V |ln(|p|/ft)r / ct 4 tt 2 



As above, for |p|/cr < l/tV/ft) 7 and (cr/ft) 1 ^ < 1, 



1 1 

< 



ln|p/ft|| " (l- 7 )|ln(a/ft)| : 
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and then, 



h = 



IMtr/fiJlV 



ajh^ 0. 



(3.25) 



Equation (3.16) follows from (3.231 and (3.25) 



We now prove (3.17). We first consider the case when |po|/f < 1- By (3.12 3.13), 



\,n. —m. I < _J_ p-<5(pL+P 2 )/2c r2 p-(Po+ 2 P-Po + (Ml-M2)Pcm-po)/o- 2 _ 1 

< ^ e- 5 (P™+P 2 )/ 2o - 2 e (Po+ 2 IPllPol+l(A'i-M2)l|Po m ||Pol)/<T 2 ( p 2 + 2|p||p | + |(Mi - ^HPcmllPoD/^ 2 < ( 3 - 26 ) 

1 e -^(pL+P 2 )/2- 2 e (l+2|p|/ CT +|Pc m |/ ff ) | po/(T | (1 + 2| p | /o . + | pcm | /(T) . 



Then, 



(ln(2+|p|/ft))° , _ s * _J_ f (ln(2+|p|/fi)) 2 ° 5(p 2 +p 2 )/(J 2 

(l + llndpl/ft)!)- 3 \Pm,p fin) 21 „4„2 J (1+| ln(|p|/7i)|)2<9 e 



e 2(l+2|p|/ ff+ |p cm |/ CT ) | po/(J |2 (1 + 2| p | /(T + \ pcm \ /a f dpcm dp . 



Estimating as in equations (3.18- 3.21) with po = 0, we prove that for |po|/c < 1, 



(ln(2 + |p|/ft))° 

(i + \H\ P \/hW 



(<An,p ~ <Pin) 



- lPn/n ° { \HaJw) ~ lP ° M ° { \H<r/K+\Po\/h)\e ) "* ' T//H - P,,|// ' lK 



In the case |po/c| > 1 the estimate is immediate from (3.15), because, 

< 



(ln(2+|p|/?»))° , 
(l + |ln(|p|/ft)|)P IWpo 



< 



IpoH o( 



|ln(<7/ft+|p Q |/R)|« 



(ln(2+|p|/ft))° 

(i+trIpTTBJIP Wpo 



as a/h + \po\/h —> 0. 



(ln(2+|p|/ft))° 

(T+TMIpTTW v?m 



(iin^/n+ipolM)!") 



□ 



We define, 



T(p 2 /2m) := S{p 2 /m) - I + in 



l + |ln(|p|/ft)l' 

where I is the identity operator on _L 2 (S X ). It follows from (2.11) and since ||<S(p 2 /2m)||g(£2(§i)) = 1, that 



Hence by (3.15 1 



|T(p 2 /2m) 



|r(p 2 /2m)v?i„, P0 || =0 



(M2+| P |/s)r 



8(L=(Si)) - w (l-(-|ln(|p|/ft)|)2' 



1 



|ln(fr/a+|po|/fi)| 



2 / ' 



asa/H+ \po\/h — > 0. 



(3.27) 



(3.28) 



(3.29) 



Let us denote, 



£(01, 02, fa, <t>i) ■■= J rfPl rfPl dp 2 dp' 2 01 (Pl, p 2 ) 02 (Pi, P2) 03(pi,P 2 ) 04 (Pi, P 2 )- 



We have that, 



7>(0) = £(0,0,0, 



(3.30) 
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The Schwarz inequality implies that, 

|£(0i, 2 , 03,^)1 <n) =1 H^-ll . (3.31) 
We state below our first result in the low-energy behaviour of the purity. 

THEOREM 3.3. Suposse that Assumption \ 1. 1\ is satisfied and that at zero H le i has neither a resonance (half-bound 
state) nor an eigenvalue. Then, uniformly for |po|/erm bounded sets, 

V(font,p )=-PCpont) + —o( 1 V asa/ft+lpol/fi^O. (3.32) 

a \\hi{o-/h+ \po\/h)\ 2 J 

Proof: Writing <p O ut, P0 as, 

</W, Po := S(p 2 /2™Vin, Po = ipin, Po - m - - * yj in , Po + T(p 2 /2m)¥V Po , 



and using (3.4), we see that we can write V(tp O \it,p ) as follows, 

4 

^(^out,p ) = 1 + £ M(P0> V'l) "02, -03, i>4) + ft(po), (3.33) 

i=l 

where 

£l,i(P0, ^1,^2,^3,^4) = £(^1,^2,^3,^4), (3.34) 

where one of the "0j is equal to T(p 2 /2TO)(^ in po and the remaining 3 are equal to <Pi n . P(r Similarly, 

A 

^•(Po) := X^^P '^ 1 '^ 2 '^ 3 '^ 4 )' ( 3 - 35 ) 
for some integer A, and where each of the £2,i(po, "0i> "02, "03j ip4) is equal to, 

C2,i(P0,1pl,1p2,ll>3,ll>4.) = £("01, "02, "03, ^4), (3-36) 

where for some 2 < k < 4, fc of the ipj are equal either to — in 1+ | EnpTSri ^".po or to 7 - (P 2 /2 m ) ( / 9 in,p and the remaining 
4 — are equal to <^i n ,p - Similarly, 

4 

Pfew) = 1 + ^i.^a.^a, ^4) + ^(0), (3.37) 

i=l 

with 

72.(0) := A(0, ^1,-02,^3, ^4). (3.38) 

i=l 

Below we prove that, 

£i,i(po, ^1,^2, Vs, V4) = £i.(0, il>i,ih,ih, ^4) + — O f M , * rTZvTa ) , i = 1,2,3, 4, (3.39) 

ft( Po ) = 11(0) + J^i O f M , /t n ,.,.,„ ) , (3.40) 

VM ff / s + Pol/a) V 



what proves the theorem in view of ( |3.33|3.37| ) 
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We proceed to prove (3.39). Without losing generality we can assume that, 



A,i(Po, ipi,ip2,ip3,ipA) = £(T(p 2 /2m)(p in , Po ,cp ilhPo ,(p in , Po ,(p intPo ). 



(3.41) 



We have that, 



£l,l(P0! tpl,tp2,1p3,1pi) = £(T(p 2 / 2m Vin, ^in,p > Vin.po > fia,po) + 
£(T(p 2 /2m)(ip in , po - Vm),^in,p 0) ¥'m,p D ,¥'m,po)' 



By (3.17 3.28 3.31 3.42 1, 



(3.42) 



Cl,l(po,1pl,1p2,1p3,1p4.) = ^(T(p 2 /2m)v5in,¥'m,po ) Vm,po,Vm,po) + ^ ° ( | i n ( ff /ft+|po|/ft)| 2 ) ' ( 3 ' 43 ) 



In the same way, using (3.17 3.29 3.43 ), we prove that, 



£l,l(P0> ^1,V>2, 1p3,i>±) = £(T(p 2 /2m)(p ia , V? in , ¥>in,p , ^m,po) + O 



a \\\n{a/h+\p \/h)\> 



(3.44) 



Repeating this argument two more times we obtain that, 



£l,l(P0)V'l>V'2,V'3,V'4) = £1,1(0,^1,^2,^3,^4) + — O ( 1 

a \\\n(a/h+\po\/h)\ 2 



(3.45) 



We prove in the same way that (3.391 holds for i = 2,3,4. Furthermore, (3.40) is proven by the same argument 



□ 



The next theorem gives us the leading order of the purity of ip ou t at low-energy. 

THEOREM 3.4. Suposse that Assumption \ 1. 1\ is satisfied and that at zero H Te \ has neither a resonance (half-bound 
state) nor an eigenvalue. Then, as a /h — > 0. 



V(<Pout) = V 



I + in £+ U7r(ln2-7+ -) - \ ) 1 

ln|p/fi| V a 2 J (In |p/^|) 2 



<Pm ) + O 



ln(a/ft)| 3 



(3.46) 



Proof: We write ip out as follows, 



^out = Vout.i + Ti(p 2 /2m)^ i: 



where, 



and 



Vout.l ; - 



I + ™ , Am S + fwr(ln2 - 7+ -) - ^ 



ln|p/ft| 



a' 2 J (In |p/^|) : 



fit 



Ti := 5(p 2 /2m) - 1 - in - ^- S- (nr(m2 - 7 + -) "" 



ln|p/fi| 



a' 2 y (ln|p//j|) 2 



By (2.11) and since ||5(p 2 /2m)|| B ( L 2( S i ) ) = 1, 



|Ti(p72m)|U,™<C 



(ln(2+|p|/ft)) 3 



B(L»(Si)) - |ln(|p|/fi)|2 



(3.47) 



(3.48) 
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Using this decomposition we write V(ip ou t) as follows, 



out.l 



(3.49) 



where TZ\{a) is given by, 



(3.50) 



for some integer D, and where each of the £j(<r, ipi,t/j2,ip3, ipi) is equal to, 



Ci(a,ipi,ip2, Ipi) = C{ipi,ip 2 , "03, -04), 



(3.51) 



where for some 1 < fc < 4, fc of the are equal to y ut,i an d the remaining 4 — fc are equal to 7i(p 2 /2m)</5; n . 
We complete the proof of the theorem proving that, 



fti(cr) = O 



ln(a//i)| 3 J ' 



as er//i — > 0. 



(3.52) 



We can assume that, 



Cx(a,1pl,1p2,1p3,1p^) = 'C(<Pout,l,^out,l,^out,l,Tl(p 2 /2m)v5in)- 



(3.53) 



By (3.16 3.481 We have that, 



£l((T, 1pl,1p2, ^3,^4) = O 



ln(a//i)| 3 



asa/H — > 0. 



(3.54) 



We estimate the remaining terms in (3.52) in the same way. 



Let us denote, 



^(q) : 



1 



e-o / 2 ,qeK 2 , 



(7T)V2 

-0in(qi,q2) := V'(qi)'0(q2) 



PROPOSITION 3.5. For any a, /? > 0, 



Proof: We follow the proof of (3.16). By (3.13) with a = 1, 



□ 



(ln|q|)« 



|in(*|qlA)l* 



o 



hy(a/h)\ 



as cr/?j — > 0. 



(3.55) 



(m|q|) Q 



< /1 + h 



(3.56) 
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where, for some 1 > 7 > 0, 
h 

and 



|q|>l/(a/fi)T 



In I 



|q|| 2a \ 



|ln( t j|q|/ft)| 2 ^ J 7T 2 



^(q c 2 m +q 2 )/ 



dq cm dq< C N (a/hy, N=l,2 



(ln|q|)° 

/|p|<l/( ff /ft)-v V I ln (cr|q|/^.)|^ y 7T 
Furthermore, for |q| < l/{a/h)'' and (oVTi) 1-7 < 1, 



1 e -^(qL+q 2 ) dqcm dq . 



< 



and then, 



|ln(a|q|/S)| " (1 - 7 )| m(a/fc)| ' 

h = O { ,, ; , tj/fi -> 0. 



|ln(a/fi)|V 



Equation (3.55) follows from (3.571 and (3.59) 



A straightforward computation with the help of (3.55) shows that, 



2 y (ln|p/fi|)2 



o 



where. 



with 



( (In(o-W )' aSCT /^0, 



Win) =^ =1 7>ljW>i„), 



and 



'Pia(V'in) = ~2tt 2 y dqidq 2 dq 3 (E^ in (qi, q 2 )) (£Vm(q3, 02)) ^m(qi) qs), 
^i : 2(^in) = -2tt 2 y dqidq 2 dq 3 (E^ in (qi, q 2 )) (£<0in(qi, qa)) ^(qa, q3), 
^1,3(^11) = 2tt 2 y dqidq 2 (S-0 in (qi,q 2 )) ^ in (qi,q2) , 

^(ipin) = 2tt 2 y dqirfq 2 (S^ in (qi, qa)) ^in(qi,q2)- 



Explicitly evaluating the integrals in (3.62 3.63 3.64) using (3.12), we prove that, 



T'ia(V'in) = — J(mi,M2), 



^1,2 



PiA^M = 27r2 (^(mi,^)) 5 

7MV>in) =2^ 2 L(/i 1 ,/i 2 ), 
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where, 



J(Mi,M2) := Jdq 2 [ /dqi Exp[-i(^ + /i^)(qi + q 2 ) 2 - (/j 2 qi - MiQ2) 2 - qf/2] 



-fodMi - a*2 I |qi + q2| lM2qi - Miq2| 

Here Jq is the modified Bessel function [31, and 



OO />00 



L(Mi,A*a):= / / dAdp e - aA e -^ + ^" (j (|Mi - wIaAp)) • 



We prove in the appendix that, 



LOi, 1 - /ii) = 



v/1 + (2 Ml - 1)2 
We denote by £(ni) the entanglement coefficient, 

£(Hi) := 2ir 2 L(p ll 1 - miti) (1 + 1 - mux)) - ^[J(pi, 1 - Ml) + JQ-- Mi>Mi)] 



(3.70) 



(3.71) 



(3.72) 



2^ 2 



1 + 71+ (2/i! - 1) 2 1 - -[J( M i,l-/ii) + J(1-mi,Mi)]- 



(3.73) 



1 + (2/ii - l) 2 

The next theorem is our main result. 

THEOREM 3.6. Suppose that Assumption \l.l\ is satisfied and that at zero H re i has neither a resonance (half-bound 
state) nor an eigenvalue. Then, 



V{<p ou t) = 1 



(Ha/h)) 



;£(»i)+0 



ln(o-/h)\ 3 J ' 



asa/H —> 0, 



(3.74) 



where the entanglement coefficient £(p,i) is given by {3.73) 



Proof: The theorem follows from $A6\ |3~60] [3~6l| [3~66pJ9l [3J2 ) 



□ 



Note that £(ni) = 5(1 — Mi)> as it should be, because V(f ou t) is invariant under the exchange of particles one and 



two. 



Observe that, 



J(1/2,1/2)=7T 3 



By (3.731 for /ii = 1/2, when the masses are equal, the entanglement coefficient is zero, 5(1/2) = 0. Of course, this 
only means that in this case the purity is one at leading order. 



We explicitly evaluate in the appendix J(1,0), 



J(1,0) = 16.6377. 



(3.75) 



For [i\ 6 [0, 1] \ {1/2, 1} we compute J(/Xi, 1 — /Ui) numerically using Gaussian quadratures. In Table 1 and in Figure 
1 we give values of £{\i\) for 0.5 < fi± := mi/(mi + m 2 ) < 1. 
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4 Conclusions 



In this paper we give a rigorous computation, with error bound, of the entanglement created in the low-energy 
scattering of two particles in two dimensions. The interaction between the particles is given by potentials that are not 
required to be spherically symmetric. Before the scattering the particles are in a pure state that is a product of two 
normalized Gaussians with the same variance a. After the collision the particles are in a outgoing asymptotic state 
that is not a product state. The measure of the entanglement created by the collision is the purity, V, of one of the 
particles in the state after the collision. Before the collision the purity is one. 

We prove that V = 1— p^^yp S + O ( | i n (a/h)\ 3 ) > as ~~ ^ ^, where a is the variance of the and the entanglement 
coefficient, £ , depends only on the masses of the particles and not on the interaction potential. This proves that the 
entanglement created at low-energy in two dimensions is universal, in the sense that it is independent on the interaction 
potential between the particles. This is strikingly different with the three dimensional case, that we considered in [3], 
where the entanglement created at low-energy is proportional to the total scattering cross section. However, the 
entanglement depends strongly in the difference of the masses. As in three dimensions |3j it takes its minimum when 
the masses are equal, and it increases rapidly with the difference of the masses. 



5 Appendix 



By H3.71D we have that [2] 

=/o D ° dAe ~ 2A i+(2 A1 2 1 -i)^ 7 ° ( i+^Vi-i? 2 ) Exp 



(2 Ml -l) 2 A 
1+(2 W -1) 2 



v/l + (2^-1)2- 



Moreover, by (3.701 and denoting q cm := qi + q 2 , 

J(1,0) = J dq 2 e~ 3q % J dq cm Exp (-q^ m + q cm ■ q 2 ) / (|qcm|M) 



Furthermore, using polar coordinates, [2], 

J(1,0) =n 3 / °° dXe- 3X 



/ °° dpe-P (lo(VpVX) 



7t 3 r o °° dXe- 2X (/ (A/2)) 2 = 7r 2 K{0.25) = 16.6377, 



(5.1) 



(5.2) 



(5.3) 



where K(x) is the complete elliptic integral. 
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Table 1: The Entanglement Coefficient £{ni) 



Hi := mij (mi + m 2 ) £{ni) 
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1.9542 
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2.1909 


0.975 


2.4216 


1 


2.6436 



5 



y 

3.0 r 




Figure 1: The entanglement coefficient y — £(p,i), as a function of x = fi\ = Wi/(roi + m^), for 0.5 < fii < 1. 
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